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Abstract 

A systematic search for superintegrable quantum Hamiltonians describing the interaction between two 
particles with spin and |, is performed. We restrict to integrals of motion that are first-order (matrix) 
polynomials in the components of linear momentum. Several such systems are found and for one non-trivial 
example we show how superintegrability leads to exact solvability: we obtain exact (nonperturbative) bound 
state energy formulas and exact expressions for the wave functions in terms of products of Laguerre and 
Jacobi polynomials. 

PACS numbers: 02.30.1k, 03.65.-w, ll.30.-j, 25.80.Dj 

1 Introduction 

The purpose of this research program is to perform a systematic study of integrability and superintegrability in 
the interaction of two particles with spin. Specifically in this article we consider a system of two nonrelativistic 
particles, one with spin s — h (e.g. a nucleon) the other with spin ,s = (e.g. a pion), moving in the three- 
dimensional Euclidean space E3. 

The Pauli-Schrodinger equation in this case will have the form 

HH> = (-^A + V (x) + ±{v 1 (x),(a,L)yjy = Ey, (1.1) 

where the V\{x) term represents the spin-orbital interaction. We use the notation 

pi = -ihd x , p 2 = -ihdy, p 3 = ~ihd z , (1.2) 

*E-mail address: wintern@crm.umontreal.ca 
^E-mail address: yurdusen@crm.umontreal.ca 
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Li = ih(zd y - yd z ), L 2 = ih(xd z - zd x ), L 3 — ih(yd x - xd y ) , (1.3) 

<ri = ( J I ) , °2 = ( ° ^ ) , ^3 = ( J \ J , (1.4) 

for the linear momentum, angular momentum and Pauli matrices, respectively. The curly bracket in (|1.1[) 
denotes an anticommutator. 

For spinless particles the Hamiltonian 

H=~A + V (x), (1.5) 

is a scalar operator, whereas H in (jl.ip is a 2 x 2 matrix operator and ^(x) is a two component spinor. 

In the spinless case (ll.5[) the Hamiltonian is integrable if there exists a pair of commuting integrals of 
motion X\, X 2 that are well-defined quantum mechanical operators, such that H, X\ and X 2 are algebraically 
independent. If further algebraically independent integrals Yj exist, the system is superintegrable. The best 
known superintegrable systems in E% are the hydrogen atom and the harmonic oscillator. Each of them is 
maximally superintegrable with 2n — 1 = 5 independent integrals, generating an o(4) and an su(3) algebra, 
respectively [TJ12]. 

A systematic search for quantum and classical superintegrable scalar potentials in (|1.5[) with integrals that 
are first- and second-order polynomials in the momenta was performed some time ago [SI El E] ■ First-order 
integrals correspond to geometrical symmetries of the potential, second-order ones are directly related to the 
separation of variables in the Schrodinger equation or Hamilton- Jacobi equation in the classical case [21 [HI H] • 

First- and second-order integrals of motion are rather easy to find for Hamiltonians of the type (|1 .5(1 in 
Euclidean space. The situation with third- and higher-order integrals is much more difficult [51 [TU1 [TTj . 

If a vector potential term, corresponding e.g. to a magnetic field is added, the problem becomes much more 
difficult and the existence of second-order integrals no longer implies the separation of variables [12 US] • 

The case (|1.1|) with a spin-orbital interaction turns out to be quite rich and rather difficult to treat system- 
atically. In a previous article we have considered the same problem in E 2 [14] . Here we concentrate on the 
Hamiltonian (jl.ip in E% but restrict to first-order integrals. Thus we search for integrals of motion of the form 

X = (A + A ■ S)p! + {B + B- a)p 2 + (C + C ■ a)p 3 + <j> + $ ■ 3 

- l Jl{(A + A-3) x + (B + B-3) y + (C + C-3) z } , (1.6) 

where Aq, Bq, Co, 4>o and Ai, Bi, Ci, <j>i (i = 1,2,3) are all scalar functions of x. 
In Section [2] we show that a spin-orbital interaction of the form 

can be induced by a gauge transformation from a purely scalar potential Vq(x) (in particular from Vq = 0). In 
Section [SI we derive and discuss the determining equations for the existence of first-order integrals. In Section 



|4]we restrict to rotationally invariant potentials Vo(r) and V\(r) and classify the integrals of motion into 0(3) 
multiplets. Solutions of the determining equations are obtained in Section [5J Superintegrable potentials are 
discussed in Section \E\ In Section [7] we solve the Pauli-Schrodinger equation for one superintegrable system 
explicitly and exactly. Finally, the conclusions and outlook are given in Section [8l 

2 Spin-orbital interaction induced by a gauge transformation 
2.1 Derivation 

In this Subsection we show that a spin-orbit term could be gauge induced from a scalar Hamiltonian (|1.5|) by a 
gauge transformation. The transformation matrix must be an element of U (2) 

e^cos(ft) e^sin(/3 3 ) \ , , 

^ -e"^ 2 sin(ft) e"^ 1 cos(ft) J ' 1 j 

where ft (j = 1, 2, 3, 4) are real functions of (x, z). It is seen that in order to generate a spin-orbit term we 
need to have 

K 2 tf(VU)- V = T (<?,£), (2.2) 

where T is an arbitrary real scalar function of (x, y, z). Equation (|2.2p implies 12 first-order partial differential 
equations for ft and T, three of which are ft^ = 0, k = 1,2,3. Hence, without loss of generality we choose 
ft = and then write the remaining 9 equations as 

cos 2 (ft)/3 M " sin 2 (ft)ft, z = , (2.3 
cos 2 (ft)ft, x -sin 2 (ftOft >:c = 0, 
fi(cos 2 (/3 3 )ft,„ -sin 2 (ft)ft., a ) - -xT, 
ft( cos(ft - ft)ft, x - i sin(ft - ft) sin(2ft)(ft + A),) - 

fr(cos(ft - ft)ft,, - i sin(ft - ft) sm(2ft)(ft. + ft)*) = xr , (2.7 

fi(sin( / 3 2 -/3 1 )/3 3 , a -|-icos(ft-/3 1 )sin(2 / 3 3 )(ft+/3 1 ) a ) = *r, (2.8 

fi(sin(/3 2 -/3 1 )/3 3 ,, + icos(/3 2 -/3 1 )sin(2/3 3 )(/3 2 +/3 1 ) z ) = -yT , (2.9 

cos(ft - ft)/?^ - i sin(ft - ft) sin(2ft)(ft + ft) a = , (2.10 

sin(ft - ft)ft :X + ~ cos(ft - ft) sin(2ft)(ft + ft), = . (2.11 

From equations (|2.6|) - (|2.11|l we obtain 

x 

Pi = V), 3 = 1, 2, 3 where £ = - , 77 = 



(2.4 
(2.5 

zl\ (2.6 



cos(ft-ft) = 3 ' £ , sin(ft - ft) = 

(ft+ft)^-^, (ft + ft), = — ^r, (2.12) 
sm(2ft) sm(2ft) 



and 



T = h }L-±- . (2.13) 

Introducing rj) into the equations (|2.3|) - (|2.5|) and using the compatibility condition for the mixed partial 
derivatives of (/?2 + we obtain the following three partial differential equations 

+ fh tVI = 2(/3 3 2 , S + Pl v ) cot(2/3 3 ) , 
/33 lf A*, - = 2(/3 3 2 4 + /? 3 2 „) (fy^ + P 3 2 ri ~ cot(2/%)/% lf ) , 

- &,„/?3,tt + /3 3)£ /?3,^ = 2(/3 3 2 c + /3 3 2 „) (r?^^ + /3 3 2 r) - cot(2&)/%,„) , (2.14) 

which could be solved for the highest-order derivatives of /3 3 (i.e. /3 3 ,fg, /3 3 ,£?j and /3 3jW ). Then, the compatibility 
conditions of these give 

^lt+Pln = - e + l 2 + v (2-15) 

which implies that r = — Hence, we conclude that Vi = A is gauge induced and it is the only potential 
which could be generated from a scalar Hamiltonian by a gauge transformation. 
The explicit form of the gauge transformation U is found as 

/3 1 = <p + c 1 , (3 2 = c 2 , f3 3 = -d + c 3 , 04 = 0, 0<6<TT, 0<ip<2n, (2.16) 

where c%, c 2 and c 3 are the following constants 

7T 

ci = C2 ± 7r , and c 3 = ±— . (2-17) 

With this transformation matrix the transformed Hamiltonian is found to be 

/ fr 2 \ fr 2 fr 2 h 

H = U- 1 i --A + Vo(x)\U = - Y A + V Q (x) + ^ + ^<j-L. (2.18) 

2.2 Integrals for V = V (r) and Vi = £ 

The potential Vi = -4- is gauge induced from a Hamiltonian of the form (|1.5p (though each term is multiplied 
by a 2 x 2 identity matrix). Hence the integrals for this case are just the gauge transforms of the integrals of 
motion of this Hamiltonian (i.e. Lj and <Xj). They can be written as 

/i fx X ' 

Ji = Li + -o-j , Si = --ci; + H-^(r, a) , (2-19) 
and satisfy the following commutation relations 

[Ji, Jj] = ihe ijk Jk , [Si,Sj] = ihe ijk Sk , [Ji, Sj] = ihe ijk Sk ■ (2.20) 
The Lie algebra C is isomorphic to a direct sum of the algebra o(3) with itself 

C ~ o(3) © o(3) = { J - S} © {S} . (2.21) 
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2.3 Integrals for V = % and V x = % 

Since these potentials are gauge induced from a free Hamiltonian, the integrals are just the gauge transforms of 
Lj , pj and <jj , which can be written as 

h h 

J l = Li + -a l , ITj = p t 2 eiklXkCTl > 




(2.22) 



They satisfy the following commutation relations 

[A - s h Sj] = o , pi, Sj] = o , [ih, iy = , 

[Ji - Si, Jj - Sj] = ihe ljk (J k - S k ) , [Ji - S t ,Uj} = ihe i:jk U k . (2.23) 

Hence the 9-dimensional Lie algebra C is isomorphic to a direct sum of the Euclidean Lie algebra e(3) with the 
algebra o(3) 

C - e(3) © o(3) = {J- S, IT} ® {S} . (2.24) 

3 Determining equations for an integral of motion 
3.1 Derivation 

In this Subsection we give the full set of determining equations obtained from the commutativity condition 
[H, X] = 0, where H is the Hamiltonian given in and X is the most general first-order integral of motion 
given in l|1.6[) . This commutator has second-, first- and zero-order terms in the momenta. By setting the 
coefficients of different powers of the momenta equal to zero in each entry of this 2x2 matrix we obtain 
the following determining equations. Since, the Planck constant h enters into the determiming equations in 
a nontrivial way we keep it throughout the whole set of determining equations. However, after giving the 
determining equations we set h = 1 for simplicity. 

i) Determining equations coming from the second-order terms 

From the diagonal elements it is immediately found that Aq, Bq and Co are linear functions and are expressed 
for any potentials Vq and V\ as 

A a = bi - a 3 y + a 2 z , B a = b 2 + a 3 x - a\z , C Q — b 3 - a 2 x + my , (3.1) 

where and bi (i = 1,2,3) are real constants. After introducing (13. 1|) into the rest of the coefficients of the 
second-order terms and separating the imaginary and real parts of the coefficients coming from the off-diagonal 
elements we are left with an overdetermined system of eighteen partial differential equations for Ai, Bi, Ci 
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(i = 1,2,3) and V\. These are, 

2zAi^i + hA 3 , x = , 2yA 1 V 1 + hA 2 . x = , 2xB 2 Vx + hB lAJ = , 2zB 2 Vx + hB 3 . y = , 
2xC 3 Vi+HC hz = 0, 2yC 3 V 1 + HC 2 , Z = , 2V X (yA 2 + zA 3 ) - = , 

2Vi (a;Si + - ftS 2 , y = Q , 2V X (xCi + yC 2 ) - W 3yZ = , 
2zVi (A 2 +Bt) + hA 3 , y + hB 3 , x = , 2yV x (A 3 + d) + HA 2 , Z + hC 2 , x = , 
2xVi (S 3 + C 2 ) + HB 1<Z + hC hy = , 2Vi (xAt + yA 2 - zCx) - hA 3 , z - hC 3 , x = , 
2V! {xBx + yB 2 - zC 2 ) - hB 3 . z - ftC 3 ,„ = , 2V X (xA 2 - yB 2 - zB 3 ) + hA ljV + %B 1>X = , 
2Vx (xAi + zA 3 - yBx) - hA 2 , y - HB 2tX = , 2V X (xA 3 - yC 2 - zC 3 ) + HA 1>Z + W hx = , 

2Vi (yB 3 - x& - zC 3 ) + hB 2 . z + hC 2 . y = . (3.2) 

ii) Determining equations coming from the first-order terms 

After introducing (|3.1j) and separating the real and imaginary parts, we have the following twelve partial differ- 
ential equations 



Vx 




- a 3 y) 


+ 2# 3 ) 


+ h(x(a Vx, x - 


V BoVx.y - 


1- Cott,,) + 


02, z , 


= 0, 


Vi 




+ a 2 z) 


- 2z0 2 ) 


+ h(x(A Vi, x - 


r B V hy - 






= 0, 


Vi 


[Hh 


- aiz) 


+ 2z4>A 


+ h(y(A V ltX - 


h B V hy - 


L C Vi,,) + 


03,xl 


= 0, 


Vi 


[h(b 2 


+ a 3 x) 


- 2x<t) 3 ) 


+ h(y{A Vi, x - 


r S F l!3/ - 


f- C0V1,,) - 


01,* 


= 0, 


Vi 


[h(bs 


- a 2 x) 


+ 2x0 2 ) 


+ k(z{A Vx, x - 




1- C 7i ja ) + 


01,J/ ; 


= 0, 


Vx 


[Hh 


+ aiy) 


- 2#i) 


+ h[z{A V^ x - 


h S Fi, a - 


h C ^ M ) - 


<t>2, x , 


= 0, 



Vx [h{a 2 y + a 3 z) - 2y4> 2 - 2z<f> 3 \ + K4>x, x = , 
Vx (h(axx + a 3 z) — 2x4>i — 2zcj} 3 J + h<fi 2 . y = , 

Vx (h(axx + a 2 y) - 2x(f> 1 - 2y0 2 ) + h^ z = , (3.3) 

<Po. x = Vx ({yA 3>x - xA 3 . y ) + (xA 2 . z - zA 2<x ) + (zAx. y - yA ljZ ) + (C 2 - S 3 )) 

+Vi,,(zA 2 - yA 3 ) + Vx, y (zB 2 - yB 3 ) + V x , z {zC 2 - yC 3 ) , 
0o,y = ^1 ((yB 3 ,i - xB 3 , y ) + (xB 2tZ - zB 2 , x ) + (zBx, y - yBx, z ) + {A 3 - C x )) 

+V ltX (xA 3 - zAx) + Vx, y {xB 3 - zBx) + Vx, z (xC 3 - zCx) , 
<f>o, z = Vx ({yC 3 , x ~ xC 3 , y ) + (xC 2tZ - zC\ x ) + {zCx, y - yC x , z ) + (Si - A 2 )) 

+V ltX (yAx - xA 2 ) + Vx, v {yBx - xB 2 ) + Vx^yCx - xC 2 ) , (3.4) 

where Aq, Bq and Co are given in (|3.1[) . There are also nine other second-order partial differential equations 

6 



for Ai, Bi, C{ and V\, coming from the coefficients of the first-order terms. However, these are differential 

consequences of (|3.2|) so we do not present them here. 

iii) Determining equations coming from the zero-order terms 

Setting the coefficients of the zero-order terms in each entry of the commutation relation equal to zero and 
separating the real and imaginary parts, we have the following four partial differential equations 

A V 0tX + B V Q ^y + CqV , z + Vi (x(4> 2 ,z ~ <f>3,y) + V{<t>3,x ~ 4>l,z) + z(4>l,y ~ <p2,x)) = , 

n{v!(Ax, z + B 2 , z - C 2 , y - Clx) + V hx ((xA hz - Mi,,) + (yA 2 , z - zA 2<y )) 
+V hy ({xB hz - zB hx ) + {yB 2 . z - zB 2 ^ y )) + V hz ({xC hz - zC hx ) + {yC 2 . z - zC 2 , y ))} 
+2(a 3 V , x + B 3 V , y + C 3 V , Z ) + 2Vx (y4>o,x - x4>o,y) = , 

(B 2 , x + C 3 ,x - A 3 , z - A 2tV ) + V lfX ((zA 3>x - xA 3 , z ) + (yA 2)X - xA 2tV )) 
+V hy ((zB 3 , x - xB 3tZ ) + {yB 2 , x - xB %v )) + V hz ((zC 3 , x - xC 3:Z ) + (yC 2>x - xC 2 , a ))} 
+2(A 1 Vo, x +B 1 V ,y+C 1 V Q , z ' S j +2V 1 (z<p ,y-y^ z ) =0, 

h^V^A^y + C 3 , y - B 3 . z - Bi.x) + Vi, x ([zA 3t y - yA 3 . z ) + (xA hv - yA hx )) 
+Vi, y ((zB 3j y - yB 3 , z ) + {xB hy - yB hx )) + Vi, z ((zC 3 , y - yC 3jZ ) + (xC hy - yC hx ))} 

+2(a 2 V , x + B 2 V , y + C 2 V , Z ) + 2V 1 [x^ z - z^ x ) = . (3.5) 

In general the partial differential equations in (|3.5|) involve second- and third-order derivatives of Ai, Bi and 
Ci, (i = 1, 2, 3), however, using (|3.2p these terms can be eliminated. 

3.2 Discussion of solution in general case 

In general the solution of the determining equations (|3.2p - (j3.5[) for the 15 unknowns <fio, Vo, V\ and Ai, Bi, Ci, 
4>i (i = 1,2,3) turns out to be a difficult problem. However, it is seen that the determining equations (|3.3p do 
not involve 4>o, Ai, Bi, Ci (i — 1,2,3) and Vo and hence could be analyzed separately. 

In order to determine the unknown functions V± and 4>i , we express the first-order derivatives of <f>i 's from (|3.3[) 
and require the compatibility of the mixed partial derivatives. This requirement gives us another 9 equations for 
</>i's and first-order derivatives of them. Now, if we introduce the first-order derivatives of 4>iS from (|3.3p into 
this system, we get a system of algebraic equations for </>,;'s {i — 1, 2, 3). This system of algebraic equations can 
be written in the following way: 

M.$ = R (3.6) 
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where M is a 9 x 3 matrix and $ and R are 3x1 and 9x1 vectors, respectively. The matrix M can be written 
as: 








2S\ 


2z5 3 









-2xd 2 


2x0$ 











2S 5 






-26i 





-2z5 4 






ZX02 


n 
u 








-2yd 2 





2yS 4 






-2z6 3 


2zSi 









~2xS 3 


-2S 6 









-2S 5 


-2yS 4 





/ 



where Si (i = 1, . . . , 6) are defined as follows 

<5i = 2Vf - 2z 2 V 2 + yV hy + xV 1}X , 5 2 = 2zV 2 + V hz , 8 3 = 2yV 2 + V hy , 

5 4 = 2xV 2 + V hx , 5 5 = 2V 1 - 2y 2 V 2 + zV hz + xV hx , S 6 = 2V 1 - 2x 2 V 2 + yV hy + zV hz . (3.8) 

The vector $ is ((f>i,(p2, <fe) T and the entries of the vector R are given as follows: 

Ru = 2V 1 (a 2 - 2xz(A V hx + B V hy + C V lyZ )) + a 2 (xV hx + yV hy + zV hz ) 

-2(b 3 x + zA )V? - b 3 V hx - z(A V hxx + B V hxy + C V hxz ) , 
R21 = y(A V hxy + B V hyy + CWuyz) + z(A V hxz + B V hyz + C V hzz ) 

-2x(h + A )V 2 + (3 - Ax^iAoV^ + B V hv + C a V hz ) - hV hx , 
R 31 = 2V 1 (a 3 + 2xy(A V hx + B V hy + C V hz )) + a 3 (xV hx + yV hy + zV hz ) 

+2(b 2 x + yA )V 2 + b 2 V hx + y(A V ltXX + B V hxy + C V hxz ) , 
R41 = -2Vi( ai + 2yz(A V hx + B V hv + C Vi,z)) - ai(xV hx + yV hy + zV hz ) 

-2(b 3 y + zB a )V 2 - b 3 V hy - z(A V hxy + B V hyy + C V hyz ) , 
i? 5 i = 2Vi(a 3 - 2xy(A V hx + B V hy + C a V hz )) + a 3 (xV hx + yV hy + zV hz ) 

-2(hy + xB n )V? - 6iVi, s - x(A V hxy + B V hyy + C a V hyz ) , 
i? 6 i = -x(A V liXX + B V hxy + C Q V liXZ ) - z(A V hxz + B a V hyz + C V hzz ) 
+2y(b 2 + B Q )V 2 - (3 - 4y 2 V 1 )(A V 1 , x + B V hy + C V hz ) + b 2 V hy , 
R n = x(A V hxx + B a V hxy + C V hxz ) + y(A V hxy + B V hyy + C V hyz ) 

-2z(b 3 + C )V 2 + (3 - 4z 2 Vi)(AoVi0 + B Vi, v + C V hz ) - b 3 V hz , 
Rsi = -2Vi(a 2 + 2xz(A Q V hx + B Q V hy + C V hz )) - a 2 {xV hx + yV hy + zV hz ) 

-2(hz + xC )V 2 - hV hz - x(AoV hxz + B V hyz + C V 1>ZZ ) , 
R gl = 2V 1 (-a 1 + 2yz(A V hx + B V hy + C a V hz )) - ai (xV hx + yV hy + zV hz ) 

+2(b 2 z + yC )V 2 + 6 2 y M + y(A V hxz + B Q V hyz + C V hzz ) , (3.9) 
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where A , B and Co are given in (|3.1| . 

In general the rank of matrix M is 3 and the rank of the extended matrix of the system (|3.6[) is 4. Thus, we 
have to analyze under which conditions we can equate these two ranks. The rank of matrix M can be 0, 2 or 3. 
It can be 2 if and only if V\ satisfies the following differential equation 

2V 1 + xV 1 , x +yV hy + zV ltZ =0, (3.10) 

which implies 

^ = ^yr^> F ^i, ( 3 - n ) 

where £ = — and rj = | . It can be if and only if V\ = rj. Hence, the rank of matrix M and the rank of the 
extended matrix can be equal if we have the folowing 

i) For rank= 0, V\ = -kt with ^ and bi ^ 0, i = 1, 2, 3 

ii) For rank= 3, V\ — V\ (r) and = for all z and a, ^ 0, i = 1, 2, 3. 
For the former case the system (|3.6p has the following solution 



11 -y 



y2 = y 



(x 2 




-^ 2 : 


)«i + 2xya 2 4 


- 2xza3 


- 26 2 ^ + 


26 3 y 








2(x 2 + y 2 - 


f z 2 ) 






(y 2 




— ar 




)a 2 + 2a;yai 4 


- 2yza 3 


+ 26iz- 


26 3 a; 








2(x 2 + y 2 - 


f z 2 ) 






(z 2 


2 

— 2T 




)a$ + 2xza\ 4 


- 2yza2 


- 2b W + 


26 2 2: 



^ 2 2(x 2 + y 2 + z 2 ) 

where ai, 0:2 and 0:3 are real constants and for the latter it has 

<Pi = y , 02 = y , 03 = y 



(3.12) 



(3.13) 



Indeed for the potential Vj. = t-, the whole set of determining equations (|3.2p - (|3.5p can be solved and we 
obtain a 9-dimensional Lie algebra C given in 1|2.24|) . For V\ — V%(r), Vq = Vo(r) we obtain the well-known 
result that H commutes with total angular momentum J = L + ^a. 

In the rest of the article we restrict to spherically symmetric potentials V\ — Vi(r), Vq — Vo(r) and have 
a rotationally invariant Hamiltonian. However, the integrals of motion would transform under the action of 
rotations 0(3). Instead of solving the whole set of determining equations we analyze the problem by classifying 
the integrals of motion into irreducible 0(3) multiplets. 

4 Classification of integrals of motion into 0(3) multiplets 

Let us assume that the Hamiltonian is 

H = -~A + V (r) + V 1 (r)c?-L, (4.1) 



and that X of (|1.6j) is an integral of motion. Rotations in E(3) leave the Hamiltonian invariant but can transform 
X into new invariants. We can hence decompose the space of integrals of motion into subspaces transforming 
under irreducible representations of 0(3). We shall also require that the subspaces have definite behavior under 
the parity operator. 

At our disposal are two vectors x and p and one pseudovector a. The integrals we are considering can involve 
at most first-order powers of p and a but arbitrary powers of x. 

General formulas for the decomposition of the representation [D(j)] n of 0(3) into irreducible components are 
given by Murnaghan [15] . Since we are interested only in j — 1 we shall proceed ab initio rather than specialize 
his results. 

We shall construct scalars, pseudo-scalars, vectors, axial vectors and symmetric two component tensors and 
pseudotensors in the space 

xpx ff\ . (4.2) 

The quantities x, p and a allow us to define 6 independent "directions" in the direct product of the Euclidean 
space and the spin one, namely 

(x, p, L = x A p, <t, a A x, a Apj , (4-3) 

and any 0(3) tensor can be expressed in terms of these. We are allowed arbitrary powers of x and any scalar in 
x space will be written as /(r) where / is an arbitrary function of r — \/ x 2 + y 2 + z 2 . Since a and p figure at 
most linearly we can form exactly 3 independent scalars and 3 pseudoscalars out of the quantities (14. 3p : 
Scalars 

S 1 = l, S 2 = (x,p), S 3 = (?,£). (4.4) 

Pseudoscalars 

P x = {a,x), P 2 = (a,p), P 3 = (x,p)(a,x). (4.5) 

The independent vectors and axial vectors are: 
Vectors 

Vi=x, V 2 =p, V 3 = (x,p)x, Vi = (a,L)x, V5 — (x, p)(a Ax), 

V 6 = aAp, V 7 = aAx, V 8 = (a,x)L. (4.6) 

Axial vectors 

Ax=L, A 2 = a, A 3 = (x,f)a, A 4 = (a,p)x, A 5 = (x,a)x, 

A & = {S,a)p, A t = (x,f)(a,£)x. (4.7) 
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Similarly we can form 10 independent 2-component symmetric tensors and 9 symmetric pseudotensors: 
Tensors 

T lk =x l x k , T 2 k = {x,p)x i x k , T^ k = (a, L) x l x k , T lk = x l p k + x V , 

if = (a,x)(x l L k +x k L l ), T lk = x\a A x) k + x k (a A x f , 

I?* = (x, p) {x l {a A x ) k + x k (a A x) 1 ) , T l & k = x l (a A p) k + x k (a A p) 1 , 

T lk = p l (a A x) k + p k (a A x) 1 , = V a k + L k a l . (4.8) 

Pseudotensors 

Yi k = (a,p)x l x k , Yi k = (a, x) x l x k , = (x,p)(a, x) x l x k , 

Y ik = (a, x) {x l p k + x k p i ) , Y* k = x l L k + x k U , Y* k = x i a k + x k a l , 

Yj k = (x,p)(x i a k +x k a i ), Y* k = p l a k + p k a l , Y^ k = L\a A x ) k + L k (a A x) 1 . (4.9) 

An arbitrary function f(r) is also a scalar and each of the quantities in (14.4p - (|4.9p can be multiplied by f(r) 
without changing its properties under rotations or reflections. 

All of the tensors and pseudotensors should be considered to be traceless since their traces appear separately 
as scalars, or pseudoscalars. For simplicity of notation we do not subtract the trace explicitly. 

5 Solution of Commutativity Equations for Vq = Vo(r), V\ = V\(r) 

In this Section we separately take the linear combinations of all the scalars, pseudo-scalars, vectors, axial- 
vectors and two component tensors and pseudotensors. For simplicity of notation we just write the bare linear 
combinations, however, in the analysis of the commutation relations we use the full symmetric form of those, 
which could be found in the Appendix. 

5.1 Scalars 

Let us take a linear combination of the scalars given in (|4.4p 

3 

It is immediately seen that in order to satisfy the commutativity equation [H, X$] = we must have 

A=ci, / 2 = 0, / 3 = c 2 , (5.2) 
where c\ and ci are real constants. The corresponding integrals <Si and S3 are trivial. 
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5.2 Pseudoscalars 

As an integral of motion we take a linear combination of the pseudoscalars given in (|4.5[) 

3 

j'=i 

and require [if, Xp] = 0. The determining equations, obtained by equating the coefficients of the second-order 
terms to zero in the commutativity equation, become 

& = -2r/ 3 Vi , h = -2/2^1 , (5.4) 
2rVf (2rVi - 3) - V/ = . (5.5) 



Depending on the solutions of the compatibility equation for V\ (|5.5p we have several cases. The solution of 
(|5.5|) is given by 

^ = i( 1 + 7T=p)- <"» 

where a 2 = 1. Note that V\ — and Vi = are special solutions of (|5.5j) with (a,/3) = (1,0) and (l,oo), 
respectively. The case V\ = tt induced by a gauge transformation has already been considered. 
Case I: Vi = -X 

For this type of potential, (|5.4I) implies that 

/ 2 = -c x r , and / 3 = — (5.7) 

r 



and the first-order equations give 



fx = - ■ (5.8) 
r 



Then zero-order equations are satisfied for any Vo(r) and we have an integral of motion Xp for the above values 
of/i. 

Casell: ^ = ^(1 + -^) , 0</3<co, a 2 = l 
For this type of potential, (|5.4I) implies that 



/2 = -^-v / l + /3r 2 , and / 3 = C ^-— , (5.9) 

P —a + v 1 + 

and the first-order equations give /1 = 0. Then Vo is determined from the zero-order equatiosns to be Vo = V\ 
and we have an integral of motion Xp for the above values of /j. 

5.3 Vectors 

Let us now take the linear combination of the vectors given in (|4.6p 

8 

Xv = Y,fo{r)V j , (5.10) 

j'=i 

12 



and require [H,Xy] — 0. The second-order terms give 



/ 2 = ci, / 3 = 0, U + fs = 0, (5.11) 

r 2 /4-/6-2r/ 4 , (5.12) 

U-2f 6 V 1 +f 8 (l-2r 2 V 1 ) = 0, (5.13) 

/£ = -2rVi(/ 8 + / 4 )-/J. (5.14) 



Equation (|5.12[) implies 
Introducing (|5.15p into (]5 . 13[) and solving for /§ we obtain 



fi= c _l±A. (5.15) 



c 2 + / 6 (l-2r 2 V 1 ) 
/8 = ^(1-2^) ■ (5 ' 16) 

Finally if we introduce / 4 and /s into (|5.14p we get a compatibility condition for V\ which is exactly same as 
(|5.5p . Thus we have the following cases: 
Case I: Vi = 7T* 

For this type of potential the commutativity equation [iJ, JCy] = implies fj = for all j — 1. ... ,8 and 
hence we do not have an integral of motion. 
Case II: Vi = ?M 1 + — | 

For this type of potential the commutativity equation [H, Xy] — implies fj = for all j = 1. ... ,8 and 
hence we do not have an integral of motion. 

Case III: Vi = Vi(r) and / 4 = , / 6 = 0, and /„ = 

For this case equations (|5.12p - (|5.14p are all satisfied and we have 

h = ci, / 3 = 0, / 4 = 0, / 5 = 0, / 6 =0, / 8 =0. (5.17) 

Then first-order terms determine /i , j-j as 

/ 1 = 0, f 7 + 2rf 7 Vx = Q, 2rf 7 V 1 + Cl V( = , (5.18) 

and we conclude that we either have V\ — ^ or we do not have any integral of motion (i.e. fj = for all 
j = l,...,8). 

5.4 Axial vectors 

Let us now take the linear combination of the axial vectors given in (|4.7p 

7 

XA = Y,M r ) A i> ( 5 - 19 ) 
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and require [H, Xa] = 0. The second-order terms give 

/i = ci, / 3 = 0, (5.20) 

/i = 4r/ 4 V r i(l-r 2 Vi), (5.21) 

/e = 2r(/ 4 - / 6 )Vi , / 4 = -/e(l - 2r 2 F x ) (5.22) 

/ 7 + 2r/ 7 Vi = , / 7 = -2/ 4 Vi . (5.23) 

Equations (f5~2T|) and (T5~2"2"f imply 

/ 6 (2Vi(3 - 6rVi + 4rV x 2 ) + rF/) = , (5.24) 

where as (|5.21j) and (|5.23[) imply 

f 4 (2rV?(3 - 2r 2 Vi) + V{) = . (5.25) 

The only common solutions of (|5.24p and (|5.25p are V\ = ■% and V\ = ^s. Thus, other than V\ = \, we have 
two cases: 

Case I: Vi = Tjpr 

For this type of potential, (|5T2Tj) - (f5"T2"3")) give 



/4 = 0, /e = J|, /7 = 0, (5.26) 



and then the first-order terms imply 



/ 2 = y, /s = 0, c 2 = 0, (5.27) 

The zero-order terms are satisfied for arbitrary Vq. Thus the only integral of motion is Xa = L + ^a, the total 
angular momentum (an integral for any V\{r) and Voir)). 
Case II: V 1 = Vi{r) and / 4 = 0, / 6 = 0, / 7 = 

For /4 = 0, fe = 0, / 7 = 0, (|5.21|) - (|5.23|) are all satisfied for arbitrary V\(r). The first-order terms give 

/ 2 = c 2 , / 5 = ( Cl - 2c 2 )Vi , (5.28) 
^ + 2r/ 5 F x = 0. (5.29) 

However, (|5.29|) together with (|5.28|) imply 

(ci - 2c 2 )(Vl + 2rV?) = . (5.30) 

Thus we have two more subcases: 
i) ci = 2c 2 

Then we have /s = and / 2 = % and the rest of the determining equations are satisfied for arbitrary V\ (r) 
and Vo(r). Hence the only integral of motion is the total angular momentum. 
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ii) V{ + 2rV 2 = 

We have V\ — r i 1 _ a ■ Introducing this V\ together with (|5.28p into the determining equations we obtain 
(ci — 2c2)a = 0. Hence we conclude that we either have a = (i.e. V\ = \) or c\ — 2c2, both of which have 
already been investigated. 

5.5 Tensors 

Let us now take the linear combination of the tensors given in (|4.8[) 

10 

*T=E-^ r ) T f> (5-31) 

and require [H, = 0. The second-order terms give 

/ 2 = 0, / 4 = 0, & = -2fi, (5.32) 

/3 + /5 + /7 + 2(/io-/ 9 )Vi=0, /3 + 2/5 + 4r 2 (/ 10 -/9)^ 2 =0, (5.33) 

h-h + 2(r 2 .h + f s + f 10 )Vi = , (5.34) 

/ 7 -2r(/ 5 -/ 7 )T4-^ = 0, (5.35) 

h + h + r 2 {h + 2(/ 10 - f 9 )Vi) = , (5.36) 

fio ~ fg + Mho - f»)Vi{l ~ r 2 V x ) = , (5.37) 

(/io - / 9 ) (2Vi(3 - Gr 2 ^ + 4r 4 Vf) + rV{) = . (5.38) 

Equation (|5.38[) implies either /io = /g or a compatibility condition for V\ which has the following solution 

Vi = — =■ ± . 1 . (5.39) 
2r 2 V4r 4 + a 

Notice that the potentials V\ = and Vi = are also solutions which correspond to limiting values of a (i.e. 
a — and a = oo). Thus we have the following cases: 

Case I: f 10 = fg 

Equations (|535]1 - (|536]1 give 

/3 = -2/ 5 , /r = /6, /s = -(/io + r 2 / 5 ), (5.40) 

for arbitarry Vi (r) and the first-order terms imply 

A = 0, / 6 = 0. (5.41) 

Then the zero-order terms are satisfied for any Vi(r) and Vq(t). However, if we introduce the above values of fi 
into the integral of motion X^? in ()5.3ip . it identically vanishes. Hence, we do not have any integral of motion 
for this case. 
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Case II: Vx = 

For this type of potential, (|5.33p - (|5.37[) imply 

fw-f 9 = j, /3 + 2/ 5 = -^, r 2 / 5 +/ 8 + /io = 0, h = .h 



However, then (15.32| implies c\ = and we are back in Case I. 
Case III: ^ = * ± * 



2r 2 \/4r' 4 +a 

For this type of potential, (|5.37[) implies 



and (|5.34[) together with (|5.36|) give 



/io = -(4r 4 + a)^+/ 9 , 



/r = ± 4cir v 1 + / 5 . 
(4r 4 + a) 4 



However, if we introduce (15.441) into (|5.35|) we see that we must have 

c\a = . 

Hence, we either have c\ = or a — 0, both of which have already been investigated. 
5.6 Pseudotensors 

Let us now take the linear combination of the pseudotensors given in (|4.9p 

i=i 

and require [if, Xy] = 0. The second-order terms give 

/3 = -2(/i+2/ 9 )Vi, / 4 --2/ 8 y 1 + (l-2r 2 y 1 )/ 9 , / 5 = 0, f 7 = f 9 , 

h = 2 (2/8^(1 - r 2 V x ) - f 9 (1 - 2r 2 y!(l - r 2 ^))) , ^ = -r(2/ 9 + r/£) , 

(/ 8 + r 2 / 9 )(2r^ 2 (3 - 2r 2 ^) - V{) +/£=(), 

(/s + r 2 / 9 )(2ry i 2 (3 - 4r 2 l/r + 2r 4 ^ 2 ) - (1 - 2r 2 V 1 )V/) = , 

(/a + r 2 h) (2rtf (3 - 4r 2 V4) + (6r 2 Vi(l - r 2 Vi) - l) V{) = . 

Equations (|5.50p and (|5.51[) imply that we either have /g + r 2 / 9 = 0, or Vi = 73- . Hence, other than V\ 
have the following case: 

Case I: / 8 + r 2 / 9 = 

Equations (|Q7|) - (|OS|) give 

/i = 2ci , h = , fi = h = h = -ci , /s = f 2 ci . 
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Then the first-order terms give 

/ 2 = 0, / 6 = 0. (5.53) 

However, upon introducing the above values of fj, (j — 1, ... ,9) into the integral of motion Ay 2 in (|5.46p . it 
vanishes identically and we do not have an integral of motion for this case. 

6 Discussion of results 

For any spherically symmetric potentials Vo(r) and V\{r) it is well known that Ji = Li + ^<7i and (&,L) are 
integrals of motion, where J is an axial vector and (<?, L) is a scalar. 

Additional first-order integrals exist only in four special cases. Two of them are treated in Sections 12.21 and 
12.31 and the spin-orbital term Vi = ^ is gauge induced. 

For Vb = V\ = -pz we obtain the vector II and the axial vectors J and S of (|2.22|) . In addition we obtain: 
Pseudoscalar 

11 i 
X P = -- (a,p) + — (a,x) (x,p) - — (a,x) . (6.1) 

Vector 

2x ._ . 2 , . . ._ . i(x A (?) . . 

V = ^ (* M + jsfr*) i^P)~ ■ (6-2) 

Axial vector 

-.1 11 x 3x 

A = —^x{o,j>) + 7f- 2 i°^)P+ ~a (°> S ) (^iP)~ i ^(^^)- ( 6 - 3 ) 

However, they all lie in the enveloping algebra of the Lie algebra (|2.24p . 

For Vi = -if and Vb = Vo(r) ^ V\ we reobtain the algebra (|2.2ip and the gauge transforms of the axial vector 
a A L and the tensor a* L k , that are in the enveloping algebra of (|2.2ip . 

For Vi = and Vq = Vo(r) we have: 
Pseudoscalars 

n=^, (6.4) 

r 

1 i 

Xp = -r(<7, p) + - (a, x) (x,p) (a, x) . (6.5) 

r r 

For Vo = Vi = ( 1 H , = I , a 2 = 1 we have the following pseudoscalar 

V yl+/3r 2 J 

Pseudoscalar 

A P = -i v ^T^( ( ?,p1 + &fi==((x,?)-i). (6.6) 
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7 An example of an exact solution of the Pauli-Schrodinger equation 

Usually the spin-orbital interaction term is treated perturbatively [16] . In the case of superintegrable systems we 
can obtain exact solutions. In this article we consider only one example: the potential Vi(r) = ot'i Vq = Vo(r) 
and the integral of motion Xp (|6.4p . The system of equations to solve is 

= , (7.1) 

J 2 * = j(j + l)*, (7.2) 

J 3 ^ = to* , (7.3) 

Xpy = e* , (7.4) 

where * = ^njmeir, 0, <f) is a two-component spinor. 
Equation (|7.3p implies: 

Equation (17. 4|) relates /i and fa: 



^MH^^AM), e 2 = l. (7.6) 
Equation (|7.2[) provides an equation for /i, namely 

h,ee + -7-/,., ~ {^(m - i)(m + | - .cos tf) - j(j + l)-\} fl = 0, (7.7) 

with 

13 5 13 . . 



2' 2 2 2 2 

To solve (|7~7| let us put 

/i = R(r)F(0) , (7.9) 

and obtain an equation for the angular part F{9) from (|7.7p . namely 

Fee + -t—^Fq - -1- (to 2 - 1 - c (to - 1) cosfl - + 1) + 1) sin 2 o\ F = . (7.10) 
sins' sin 9 ^ 4 2 V 4/ J 

Equation ([7T0|) can be solved in terms of Jacobi polynomials. We get different expressions for to > and to < 0. 
For to < we have 



F(9) = (l-z 2 )i-^P^\z), (7.11) 

OL = — TO -)- ^ , /3 — — TO — — j 2 = 0080, (7-12) 
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which is regular for — 1 < z < 1 (for m < 0). 
For m > we have 

F(0) - (1 - z)*-4 + i(l + + Pj"^^) , (7.13) 

Of = 771 — — , /3 = 777 + - , 2 = cos# . (7-14) 

Solution (|7.13p is regular for — 1 < z < 1 (for m > 0). 

Finally, to obtain the radial part of the solution we put the results obtained so far into (|7.1[) and obtain the 
radial equation: 

\ { R " + l R ) + { Vo{r) + ( J{J + l)^s} R = ER - (7 ' 15) 

We shall solve (|7.15[) for the case when Vo(r) is the Coulomb potential 

Vo(r) = -, M<0. (7.16) 

r 

The result is obtained in terms of Laguerre polynomials. We put 

R(r) =e wr r p L{ar), (7.17) 

and obtain 

orL" + 2(p + 1 + -ar)£' + (— fp(p + 1) - j(j + 1) + + 2(u>(p + 1) - /j) + — ^^ arl L = . (7.18) 

a yor \ A J a a ) 

This coincides with the equation for Laguerre polynomials L = L" (or) if we put 



1 /.„ 1 



-(«>(p+l)-/i) =», OL=2sjp+j--. (7.19) 
From (|7.18p we also find the bound state energies to be 

£ nj = ; ^ , , . (7.20) 



We see that the energy depends on only two quantum numbers, n and j whereas the wave function (|7.5p 
depends on n, j, m and e. The spin-orbital interaction removes the "dynamical" or "accidental" degeneracy 
with respect to the quantum number j. Superintegrability relates the two components of the spinor $„j me and 
this made it possible to calculate the wave functions explicitly and exactly. 
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8 Conclusions and outlook 



The main results of this article can be summed up as a theorem. 

Theorem 1. The only first-order spherically symmetric superintegrable systems of type Q1.1JI in the Euclidean 
space E3 are the following ones: 

ti 2 h 

1. % = Vi = £. (8.1) 

The integrals of motion are given in (|2.22[1 and form the algebra (|2.23|) , (|2.24p . The potentials (|8.ip are induced 
from free motion by a gauge transformation. 

2. V = V (r), Vi = ^, (8.2) 

where Vo(r) is arbitrary. The integrals are given in (|2. 19|) and form the algebra (|2.20j) . 12.21 j) . The spin-orbital 
term V\ is induced by a gauge transformation. 

3. V = V (r), V i = ^2, ( 8 - 3 ) 
where Vo(r) is arbitrary. The integrals are the two pseudoscalars (|6.4[) and (|6.5[) . 

4. V = hV u y 1 = A(i+ a V a 2 = i. (8 .4) 

The integral is the pseudoscalar (|6.6|) . 

In all cases the components of the total angular momentum Jj = X, + |o"i, i = 1,2,3 are also integrals of 
motion. □ 

In the case of scalar particles (Vi(r) = 0) first-order superintegrability does not exist. The best known 
cases of second-order superintegrability are the Coulomb atom and the harmonic oscillator. In the first case the 
additional (to angular momentum) integrals of motion form a vector (the Laplace- Runge-Lenz vector). In the 
second case they form a two-valent tensor. In both cases the integrals generate a non Abelian Lie algebra and 
this leads to an additional degeneracy of the energy levels. 

In the case considered in this article the situation is different. First of all, first-order superintegrability does 
exist (see Theorem 1). For cases 3 and 4 of the Theorem the additional integrals are pseudoscalars and hence 
commute with the total angular momentum. Hence it is possible to simultaneously diagonalize H, J 2 , J3 and 
the additional integral X. In the example, considered in Section[71 the energy depends on two quantum numbers 
n and j, the wave functions on n, j, m and e = ±1. We thus have the geometric degeneracy related to the 
operator J3 and also a discrete degeneracy due to X . 

Theorem 1 also provides examples of "pure quantum integrability" [TU1 HH E]- The potentials V\ and 
sometimes also the potentials Vq disappear in the classical limit h — > 0. 
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It has been conjectured [18] that all maximally superintegrable (scalar) systems are also exactly solvable. 
This means that their bound state energies can be calculated algebraically. Moreover their wave functions can be 
expressed as polynomials in the appropriate variables, multiplied by some overall factor. The example (|7.ip - (17.4|) 
is superintegrable, but not maximally so. We have however shown that for Vq = ^ the system is exactly solvable. 

The conjecture of Ref. [T5] has been supported by many examples [TBI ITOI |2T)] . 

In a future article we plan to study the potentials (|8.3[) and (|8.4j) in more detail, making other choices for Vq 
in (|8.3p and diagonalizing a more general operator X2 + aX\ (see (|6.4[) and ()6.5ll ). 

Another project that is being pursued is that of second-order superintegrability. The Hamiltonian is the 
same as in (ll.ip . however, the integrals, additional to total angular momentum are not of the form (|1.6[) but are 
second-order polynomials in the linear momentum p. 
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APPENDIX 

In this Appendix we give the full symmetric form of the integral of motions separately for scalars, pseudoscalars, 
vectors, axial vectors, tensors and pseudotensors. 

i) Scalars 

The full symmetric form of X$ is given as 

X s = fi + h (x,f) + h Z)-\r&-%h- (A-l) 

ii) Pseudoscalars 

The full symmetric form of Xp could be given as 

X P = h (<?, x) + h (3,0) + h [3, x) (x,p)-i ^ [f' 2 + r 2 /a + 4r / 3 ) . (A-2) 

iii) Vectors 

The full symmetric form of Xy can be written as 

Xv = a ? (/i-^(y+r^ + 4/ 3 )+/ 3 (2 ? ,p) + /4(a,i)) 
+ f2P + f6(3Ap)+f s (a,x)L 

- ^Af)(^+/4-/ 8 +^+2/ 5 + 2i(/ 5 (i,p)+/ 7 )). (A-3) 

iv) Axial vectors 
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Let us now give the full symmetric form of Xa 

Xa = fiL + a (/ 2 - l - (3/ 3 + rf 3 + h + k) + h % P)) 

+ x (a, x) (/ B - £ (fi + f 6 ) l - (5/ 7 + r/ 7 ) + / 7 (af, p)) 

+ f 4 x(a,p) + f 6 (a,x)p. (A-4) 

v) Tensors 

In the commutator relation [H, \ — 0, it is enough to consider only one component since the others then 
necessarily commute due to the rotations. Let us now give the full symmetric form of X)j? 

X? = *y(h + f2(x,p) - l -{rf> 2 + 5/2) + fs{ff, L) - i&) 

+ (zxat - zya 2 - (x 2 - y 2 )a 3 ) (1 (/ 3 - f 5 + rfi + 5/ 7 + \ (/£ + &)) - f 6 - f 7 (x, p)) 

+ U{xp y + yp x ) + f 5 (a, x)(yL 1 + xL 2 ) + h{y<J2Pz ~ yo- 3 p y + xa 3 p x - xa x p z ) 

- h ((y&3 - z<T 2 )p y + (zai - xa 3 )p x ) + fi (<T 2 L 1 + 0-1X2) • (A-5) 

vi) Pseudotensors 

Similarly it is enough to consider only one component. The full symmetric form of Xy is 

X Y 2 = hxy(a,p) - xy(a, x) (*£ _ f 2 _ frfafl + % r A + 3i f 3 + 

- (yai + xo 2 ) + U + rfr - h) - fe - fr&ft + 2iA + i^) 

+ U(a, x)(xp y +yp x ) + h{yLi + xL 2 ) + fo(<J\p v + o 2 p x ) 

+ f 9 ([xa 3 - zai)Lx + (za 2 - y(T 3 )L 2 ^ ■ (A-6) 
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